High-Precision Limb-Darkening Measurement of a K3 Giant Using 

Microlensing 

Dale L. Fields 1 
and 

M. D. Albrow 2 * 3 , J. An 1 ' 4 , J.-P. Beaulieu 5 , J. A. R. Caldwell 6 , D. L. DePoy 1 , 
M. Dominik 7 , B. S. Gaudi 8 ' 9 , A. Gould 1 , J. Greenhill 10 , K. Hill 10 , U. G. J0rgensen n , 
S. Kane 7 , R. Martin 12 , J. Menzies 6 , R. W. Pogge 1 , K. R. Pollard 2 , P. D. Sackett 13 , 
K. C. Sahu 3 , P. Vermaak 6 , R. Watson 10 , and A. Williams 12 

(The PLANET Collaboration) 

and 

J.-F. Glicenstein 14 , and P. H. Hauschildt 15 
ABSTRACT 

We obtain high-precision limb-darkening measurements in five bands (V, Ve, Ie, I, 
and H) for the K3 III (T eff = 4200 K, [Fe/H]= +0.3, log g = 2.3) source of the Galactic 
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bulge microlensing event EROS BLG-2000-5. These measurements are inconsistent with 
the predictions of atmospheric models at > 10 a. While the disagreement is present in 
all bands, it is most apparent in /, Ie and Ve, in part because the data are better and in 
part because the intrinsic disagreement is stronger. We find that when limb-darkening 
profiles are normalized to have unit total flux, the /-band models for a broad range 
of temperatures all cross each other at a common point. The solar profile also passes 
through this point. However, the profile as measured by microlensing does not. We 
conjecture that the models have incorporated some aspect of solar physics that is not 
shared by giant atmospheres. 



1. Introduction 

The brightness profiles (limb darkenings) of stars are a potentially powerful probe of their 
atmospheres as a function of depth. At each point along the projected radius of a star, the observed 
flux originates from a range of physical depths, the deepest of which (the surface of last scattering) 
increases with physical radius as one's line of sight progresses from the center of the star towards 
its limb. Hence, since stellar temperatures generally fall towards the surface, one expects that the 
limb will appear cooler (and therefore redder and fainter) than the center. If model atmospheres 
accurately reflect the physical conditions of the star as a function of depth, they should reproduce 
the star's limb-darkening profile. 

Because limb darkening is a photometric quantity, it can, in principle, be measured to high 
precision. The drawback is that one has to be able to determine where on the star the light is 
coming from. Historically, there are two ways this is done. The first is resolving the star. The most 
obvious example would be the Sun (Pierce & Waddell 1961). Recent advances in interferometry 
have allowed one to resolve the surfaces of the highest angular-diameter stars (Burns et al. 1997) 
and have in recent years provided data good enough to begin challenging models. The second 
method is by occultation, either by an object in our solar system or one orbiting the observed star. 
The Moon is the only occulter used in our solar system. However, while lunar occultations are 
sufficiently precise to demonstrate that limb-darkened models are superior to uniform-brightness 
models, they lack the precision to test limb-darkening models. If the occulting body is in the source's 
system, it can be a star or a planet. If it is a star, the system is more properly referred to as an 
eclipsing binary. While such systems would seem to have great potential, it is extremely difficult 
to disentangle the limb-darkened profile from other parameters describing the fit to an eclipsing 
binary lightcurve (Popper 1984, 1985). The first extra-solar transiting planet to be discovered is 
HD 209458b (Charbonneau et al. 2000). Because the planet is much smaller and darker than the 
star, its transits can be used to trace the stellar light profile in great detail (Brown et al. 2001). 
Moreover, it is generally expected that ongoing and future transit surveys will turn up several more 
such systems. 
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There is one other method that can distinguish between light coming from different parts of a 
star. If a star passes through a microlensing pattern, different parts of the star are magnified by 
different amounts. In practice, differential magnification is significant only when the star passes 
through a caustic, which is a region of formally infinite magnification for a point source. Up to now, 
the best measurement of microlens limb darkening has been from MACHO 97-BLG-28 (Albrow et 
al. 1999a). This is because the event included a cusp crossing and thus the magnification pattern 
was sharper, thereby giving better resolution across the star. In addition, the data for this event 
are very good. Albrow et al. (1999a) were able to measure two limb-darkening parameters each 
in V and / for the K giant source. They demonstrated that the resulting surface profiles are in 
reasonable agreement with the predictions of atmospheric models for stars of the same spectral 
type. However, they were not able to challenge these models. 

Afonso et al. (2000) obtained a linear limb-darkening coefficient in each of four bands for an A 
star in the Small Magellanic Cloud (MACHO 98-SMC-l). They confirmed the expected trend of 
increased limb darkening toward the blue, but the measurements are not precise enough to challenge 
models. The primary difficulty is that the source star is extremely faint, I ~ 22, so that even when 
it is highly magnified, the signal-to-noise ratio (S/N) is modest. 

Albrow et al. (2000a) obtained linear limb-darkening coefficients for a red clump giant (MA- 
CHO 97-BLG-41). Even though the event itself was quite favorable, with three caustic crossings 
and a cusp crossing, bad weather and bad luck combined to limit the sensitivity of the data to limb 
darkening. 

The limb darkening analysis of OGLE-1999-BUL-23 by Albrow et al. (2001a) was a major 
breakthrough in this subject. They developed a method to simultaneously compare limb-darkening 
measurements in two bands (V and /) to the predictions of a whole suite of atmospheric models. 
The analysis demonstrated a conflict only at the 2 a level, so no significant conclusions could 
be drawn. However, if the density of measurements had been higher or the errors smaller, this 
technique would have been able to give observational input into atmospheric modeling of the limb 
darkening of a moderately evolved star for the first time. 

The binary-lens microlensing event, EROS BLG-2000-5 provides the best constraints on limb 
darkening by any microlensing event. This event has a caustic crossing that is four days long. 
This extraordinarily long timescale and the generally excellent weather for all four days at all 
five observatories combine to yield an extremely high density of coverage of the source crossing 
in units of its own radius. The first caustic crossing is well measured (fortuitously since its onset 
cannot be predicted beforehand), and the event contains a cusp approach in addition to the two 
caustics. Hence, this event is better constrained than any other microlensing event (An et al. 
2002). With the physical parameters of the event well constrained, higher order terms in the 
microlensing parameterization (such as limb darkening) can be precisely determined. Finally, we 
have spectroscopic data for this star (taken when it was highly magnified) that give us independent 
information about its temperature, metallicity and surface gravity. This combination of information 



-4 - 



allows us to to confront model atmospheres of what we determine to be a K3 III star in all five 
bands. We will test the Kurucz ATLAS models (Claret 2000) in the Johnson-Cousins V, I and H 
bands, and the Hauschildt Next 2 Gen models in the same Johnson-Cousins bands plus two EROS 
bands, Ve and Ie- The Next 2 Gen models are the current versions of the NextGen models described 
in Hauschildt et al. (1999a,b). 

2. Data 

In our analysis of EROS BLG-2000-5 we make use of 11 data sets in five filters. The PLANET 
collaboration contributes nine data sets in three standard filters, 3 sets in V, 4 in / and 2 in H. A 
description of the / band data can be found in An et al. (2002). The V band data are very similar 
in quality to the I band, the main difference being that they contain about half the number of 
points. The H band data were taken at SAAO and YALO by the instruments DANDICAM and 
ANDICAM, respectively. The instruments and procedures are identical in these two cases, and 
each contains a Tek 2048 x 2048 CCD and a Rockwell 1024 x 1024 HgCdTe IR array. The light 
path contains a dichroic that allows optical and near-IR images to be obtained simultaneously at 
the same position on the sky. The H band images are constructed by averaging five contiguous 
dithered frames of 60s each. The dithered images are flat-field corrected and then used to create 
a median sky image, which is subtracted from the individual frames before they are shifted and 
co-added. The last two data sets are from the EROS collaboration. These consist of V-like and 
/-like bands (Ve and Ie) that are described in Afonso et al. (2001). In the same manner as An et 
al. (2002), each data set is cleaned of bad points and has its errors rescaled so that the reduced x 2 is 
unity at the best-fit model. Note that the number of points in the PLANET / band data sets may 
be individually different from that given in An et al. (2002) as these cuts were done independently. 
However, since the total difference is only 10 points out of 1287, the impact on our conclusions is 
negligible. The attributes of the data sets are given in Table 1. 

3. Model 

We continue the model formalism of An et al. (2002), which contains 11 geometric parameters. 
Seven of these are static binary lens parameters: the lens separation in units of the Einstein radius 
d tc (henceforth simply d), the binary lens mass ratio q, the angle between the direction of motion 
of the source and the binary lens axis a' , the distance between the cusp and the source at closest 
approach u c , the time taken to travel an Einstein radius t' E , the time of closest approach to the 
cusp t c , the ratio of source radius to Einstein radius p*. Two are rotational parameters: d and 
ui, and two are vector components of microlens parallax: 7Te,|| and tte,±- The derivation of this 
parameterization and its relation to the standard formalism is given in An et al. (2002). Our 
geometrical solution and that of An et al. (2002) are given in Table 2. 
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In addition, each observatory and band has its own five photometric parameters: the unmagni- 
fied source flux f s , the blend flux fa, a linear seeing correction rj s , a linear limb-darkening parameter 
r, and a square-root limb-darkening parameter A. These photometric parameters are returned from 
a linear fit to the magnification curve determined by the 11 geometrical parameters. For each band 
(V, I, and H) that is observed from several observatories, all observatories are constrained to give 
the same values of V and A. 



The form of the limb-darkening law we use is, 



Sx(#) = S x 



(1 — T\ — A A ) + — 7J— cos i? + ^ V^d 



(1) 



which conserves flux independent of T and A, with S\ being the mean surface brightness of the 
source and •& the angle between the normal to the stellar surface and the line of sight. This law is 
a different form of the more widely-used, 



S x (tf) = S x (0) 1 - c A (l - cos i?) - d x (l - Vcostf) . (2) 

It should be noted, however, that equation (2) is normalized to the flux at the center, and thus the 
total flux is a function of S\(0), c\ and d\. 

The transformation of the coefficients in equation (1) to the usual coefficients used in equa- 
tion (2) is given by, 

- 6r A , _ 5A A . . 

CA 4 + 2r A + A A ' ax 4 + 2r A + A A {6) 

while the inverse transformation is r A = 10c A /(15 — 5c A — 3d A ) and A A = 12d A /(15 — 5c A — 3d\). The 
limb-darkening parameters are primarily determined by the behavior of the light curve between the 
time the source edge enters the caustic and the time the source center enters the caustic (and the 
inverse of this process as the source leaves the caustic). 

We then take this model and expand it to include all five bands, constraining all observations 
in the same band to give the same limb-darkening parameters. However, since there are no seeing 
data for EROS, these two bands do not have seeing corrections. This gives our model a total of 11 
(geometric) + [5 (photometric) x 11 (data sets) —2 (EROS seeing coefficients)] = 64 fit parameters, 
which are then subject to [(2 V + 3 I + 1 H ) x 2 (limb-darkening parameters)]= 12 constraints. 



4. Analysis 

4.1. x 2 Minimization 

We then take this model and search for minima in \ 2 i n t ne 1 1-dimensional space of the geo- 
metrical parameters. We employ a simple grid search algorithm and not a more efficient technique 
such as simplex because the (apparent) \ 2 surface is rough with many false minima. This problem 
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forces us to restrict our automated grid search to nine of the geometrical parameters and to step 
through the remaining two (d,q) "by hand." 

Even so, we find that we can locate the true minimum only to within ~ 0.02 in d and q, 
despite the fact that the true errors in these parameters (as determined from the curvature of the 
X 2 surface measured over larger scales) are < 0.01. That is, the apparent \ 2 varies by ~ 10 for 
the same (d,q) when we initialize our search using different values of the other nine parameters. 
We believe that this roughness is most likely due to numerical noise, rather than roughness in the 
"true" x 2 surface. However, regardless of the exact cause of the roughness of the surface, its impact 
via increased uncertainty in (d,q) on the errors in the limb-darkening parameters must be assessed. 
This will be done in §4.2. 

We start the grid search near the minimum found by An et al. (2002). To save computation 
cycles, we first do a grid search using only /, and H band data. After the new minimum has 
been approximately located, the V and Ve band data are finally used. 

As stated above, the fitting routine finds both linear T and square-root A limb-darkening 
parameters. However, in the case of the H band, the degeneracy between Th and Ah is too 
severe for this 2-parameter fit to provide useful information. Because Th is consistent with zero, 
we constrain it to be exactly zero and report only Ah- The analysis is then redone giving us 64 
parameters subject to 13 constraints, which, since two parameters are constrained to be zero, is 
equivalent to 62 parameters with 11 constraints. The limb darkening of the best-fit microlens model 
is shown in Figure 1. 

4.2. Errors 

Contributions to the error in the limb-darkening parameters can be broken down into three 
sources: the photometric (which include limb-darkening) parameters, the geometric parameters we 
minimize over, and d and q. The covariance matrix from the photometric parameters is easily 
obtained as it is a byproduct of the linear fit that solves for the photometric parameters. This is 
done holding all the geometric parameters fixed. We then apply the hybrid statistical error analysis 
as given in Appendix D of An et al. (2002) to determine the combined covariance matrix for the nine 
geometrical + 53 photometric parameters. While this approach does yield an invertible covariance 
matrix, there are reasons not to trust this method in this particular instance. We have stated above 
in §4.1 that the apparent x 2 surface is rough, with numerical noise a possible culprit. An effect of 
numerical noise upon the ideal x 2 parabola is to artificially raise the x 2 at any particular point in 
parameter space. This is illustrated in Figure 2. A good approximation to the true (numerical-noise 
free) parabola can be found by fitting to the envelope as noise will not decrease x 2 - We fit for the 
error induced by the nine geometric parameters that we minimize over in this fashion. We fit for 
the error in each band separately. In all cases, we find that the T-A error ellipse induced by the 
nine geometrical parameters is small compared to that induced by the photometric parameters. 
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The two sets of error ellipses are well- aligned, differing in direction by only a few degrees. We add 
the resulting covariance matrices to obtain the geometric + photometric errors. 

Up to this point, we have not yet taken into account the error induced by d and q. The 
roughness of the x 2 surface is also a factor in this analysis. Unfortunately, while it is theoretically 
possible to use the same method that we use on the other geometrical parameters, it would take 
several orders of magnitude more computational resources than we currently have to properly 
populate the geometrical figure analogous to Figure 2. We instead investigate whether and to what 
extent the additional error in our limb-darkening parameters induced by d and q will affect our 
conclusions. 

To this end, we first define 

(A x 2 )LD = ^2Sa i b ij 5a j , (4) 

i,j 

where 5ai is a vector whose components are the differences between the limb-darkening parameters 
of our best-fit model and a comparison model. Here, bij = (cjj) -1 and Cij is the covariance matrix 
of the limb-darkening parameters evaluated at the best-fit model. We define equation (4) as A% 2 
because it is the distance, expressed in x 2 and normalized by the error ellipsoid, between models 
that both attempt to describe the data. We then compare our best-fit microlens model with other 
microlens models over the roughly (0.02 x 0.02) region of the (d,q) space over which \ 2 cannot be 
properly minimized on account of the roughness in the x 2 surface described in §4.1. We find that 
the limb-darkening parameters vary by less than 1 a. By comparison, as we will show in §4.4 and 
§4.5, the limb-darkening parameters of our best-fit model differ from those predicted by the stellar 
models by > 10 a. We conclude that the additional errors resulting from both the uncertainty in d 
and q and the underlying numerical noise do not affect our overall conclusions. We recognize this 
additional error exists, but given that we have no way to quantify it properly, we simply report the 
error induced by the photometric and other nine geometric parameters in Table 3. 



4.3. Independent Analysis of Source Star 

When we compare our limb-darkening results with the atmospheric models in §4.4 and §4.5, 
we wish to restrict attention to models that are relevant to the source star. We therefore begin 
by summarizing the results of an analysis of the source's physical properties as given by An et al. 
(2002). Assuming no differential reddening across the field, the dereddened color and magnitude 
of the source can be found by measuring the source offset from the red clump identified in a color- 
magnitude diagram (CMD) of the field: (V-I) = 1.390 ± 0.010 and I = 14.70 ± 0.03. From its 
color and the fact that it is a giant (see below), the source is a K3 III star with a corresponding 
T e ff = 4200 K. The color and magnitude imply a source angular radius 0* = 6.62 ± 0.58 //as. The 
source-lens relative proper motion and source radial velocity imply that the source most probably 
lies in the bulge, i.e., at d s ~ 8kpc, which implies a physical radius r* = d s #* = 11.4 ± 1.0 R Q . 
From kinematic information, An et al. (2002) find that while it is possible that the source lies in 
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the far disk (and so is physically bigger), it essentially cannot lie in the near disk. Since the source 
is either in the bulge, or in the far disk ~ 500 pc above the plane, it must be a fairly old giant and 
therefore have a mass M ~ 1M Q . The mass and radius combine to give log g ~ 2.3. If it lies in 
the far disk at 12kpc then log g can be as much as 0.3 dex smaller. 

Ramirez et al. (2003) have analyzed HIRES Keck spectra taken on the last two nights of the 
caustic crossing, HJD 2451731.953 and HJD 2451732.950 when the source center was approximately 
0.25 and 0.75 source radii outside the caustic. By comparing these spectra with integrated spectra of 
model sources, they estimate that T ef j = 4250 K, \ogg = 1.75, and [Fe/H]= +0.29 + 0.04 on the first 
night and that T cS = 4450 K, log g = 2.25, and [Fe/H]= +0.22 ±0.07 on the second. The gravity is 
too loosely constrained in this analysis to be of any use, but fortunately we have the photometric 
method applied to the An et al. (2002) results given above, which is both simpler and more robust. 
Because the source is differentially magnified while the models are not, this spectroscopic approach 
is not fully self-consistent. Nevertheless, we expect the error induced to be modest, particularly 
on the first night when the limb of the star is not particularly emphasized in the integrated source 
light. Since, in addition, both the observing conditions and S/N were substantially better the first 
night, we adopt T e g = 4250 K and [Fe/H]= +0.3 as the spectroscopic determinations. The former 
is in excellent agreement with the photometric determinations described above. We designate 
(T eff ,[Fe/H],log#)=(4200K,+0.3,2.3) as the most physical model (MPM). 

If we are to use these estimates to define a viable region of model atmosphere parameter 
space, we need error estimates as well. Spectroscopic temperature estimates are routinely good 
to a hundred Kelvins. Similarly, our photometric temperature estimate can be off by a hundred 
Kelvins depending on differential reddening. As summarized above, the source gravity is strongly 
constrained by the angular-size measurement and distance estimates. One could possibly push the 
source into the near part of the bulge, or into the far disk, but that is all. There is also some 
error associated with the mass estimation, but due to age constraints the mass cannot be too far 
from ~ IMq. To be safe, we budget a 50% mass error. Such errors give us ±0.38 in log g. The 
metallicity is the least well constrained. This is only a minor problem, for as we shall show in §4.4 
and §4.5, metallicity has only a minor effect upon the model-atmosphere limb-darkening curves. 
We set our lower metallicity limit at solar. Physically reasonable models (PRMs) would then be 
T eff = 4100-4300 K, logs = 1-9-2.7 and [Fe/H]=0.0 - +0.3. 

4.4. Comparison with ATLAS Models 

Claret (2000) fits five different limb-darkening laws to a suite of ATLAS model atmospheres 
supplied to him by R. Kurucz in 2000. Claret (2000) then reports the parameters for each of these 
limb-darkening laws. We compare our results to the linear + square root law, rather than the favored 
four-parameter fit since we also use a linear + square root law in our fitting and thus the coefficients 
are comparable. Little is lost by this substitution as the four-parameter and two-parameter fits 
differ by much less than the difference between the microlensing-based and atmospheric-model 
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profiles. We use equation (4) to create our measure of goodness of fit, but because we are now 
comparing theoretical atmospheric models to what we consider a parameterization of reality, we 
term the result x 2 ■ Eq. 4 implicitly assumes that the x 2 surface is parabolic which, since the 
microlensing fit is non-linear, is not strictly the case. However, as shown in §4.2, the covariance 
matrix is dominated by the linear part of the fit. Hence the x 2 surface is nearly parabolic. As 
before 5a>i is a difference in the limb-darkening parameters, this time between an atmospheric 
model and our microlensing result, and bij is the inverse of c^, the covariance matrix for the 
microlensing limb-darkening parameters. We convert the Claret (2000) c\ and d\ to T\ and A\ 
using the inverse of equation (3) to make the comparison. Since we have enforced Tjj = on the 
microlensing models, we repeat this constraint on the atmosphere models. To do this, we determine 
the (r#,A/f) covariance matrix by removing the constraint Th = in the microlensing fit. We 
then use the method given in Appendix A of Gould & An (2002) to enforce the constraint on the 
Claret (2000) data. We restrict the comparison to the standard (Johnson-Cousins) bands, V, I, 
and H, because Claret (2000) reports limb-darkening parameters only for these. 

Following the procedure pioneered by Albrow et al. (2001a), we begin by simultaneously com- 
paring the microlensing limb-darkening parameters from all three filters with the atmospheric-model 
parameters, taking full account of the covariances among these five parameters. We restrict our 
investigation to those models with turbulent velocity vt = 2 only. Models given by Claret (2000) 
with other vt do not span the full parameter space required by our investigation, nor do we have 
independent information that would distinguish among different vt as we do for T e g, [Fe/H] and 
log g. Special attention is paid to two regions of the x 2 surface: the neighborhood of the MPM to 
check for consistency between the atmospheric models and the microlens data, and features around 
Xmin (which may not be near the MPM) to try to guide modelers in understanding the results 
of their simulations. The ATLAS parameter grid does not contain the MPM, but the closest is 
T eff = 4250 K, log g = 2.5 and [M/H]= +0.3. We shall refer to this as the MPM while within §4.4. 
Additionally, the ATLAS models given in Claret (2000) have a larger grid spacing than the region 
covered by the PRMs, so we investigate T cS = 4000-4500 K, log g = 2.0 and 2.5, and [M/H] from 
0.0 to +0.3. We define consistency as having a x 2 ^ 4. We find that x 2 a t the MPM is 164. In fact, 
all of the PRMs are high, with the x 2 minimum at 45. This in itself is a major concern. That no 
model atmosphere agrees with our data, regardless of its parameters is evident by the fact that the 
x l iin = 37. This occurs at T cS = 4500 K, logg = 3.0 and 3.5, and [M/H]= -0.3, which is incom- 
patible with the other evidence we have about this star. We find that, in all the ATLAS models 
in the vicinity of the MPM, the differences in x 2 between log g = 2.0 and 2.5 are small compared 
to those induced by changes in the other two parameters. This is true for the / and H bands as 
well, but log g has more of an effect in the V band. Due to this lack of distinguishability, we shall 
focus our investigation on log g = 2.5. Any model whose gravity is not listed should be assumed to 
have log ^ = 2.5. To determine whether the large mismatch in limb-darkening parameters comes 
primarily from one specific band, we investigate the goodness of fit for each band separately. 
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44-1- V band 

We repeat our x 2 minimization over the space of ATLAS models considering only the V band 
parameters Ty and Ay. The MPM has a x 2 of 19. We then look for PRMs that might be consistent. 
All points with T c g = 4500 K have x 2 < 4, and are thus consistent. Within that consistency, higher 
[M/H] and log g are favored. Remember, however, that the grid spacing of the ATLAS models is 
larger than the true permitted temperature range. If we take this into account, and note that at 
T e s = 4250 K, x 2 > 11) we must downgrade V band to marginal inconsistency. The shape of the x 2 
surface near Xmin re turns a "valley" (part of which is shown in Fig. 3) running along T c g = 4500 K. 
Surface gravity is more important in the V-band x 2 than in the sum over all the bands. Regardless, 
X 2 in the V band alone still depends more upon [M/H] than log g. 

44.2. H band 

We perform an analysis for H band in the same manner as in the previous section. The MPM 
has x 2 = 16, similar to that of the V band. The H band also has consistency at the same level and 
manner as V band, at T e g- = 4500 K in the PRMs and the same caveat discussed in §4.4.1 applies 
for the most part. The H band also has a valley structure (Fig. 3) in its x 2 surface analogous 
to that in the V band, though there is some dependence on T c g in the range that we investigate. 
This T e s dependence is slight; the location of the "valley floor" only shifts from T e g = 4500 K at 
[M/H]= -0.3 to T cS = 4750 K at [M/H]= +0.3. As in the all-bands, \ogg is essentially unimportant 
in the H band. 

44.3. I band 

The I band MPM has x 2 = 128. None of its PRMs are consistent with the microlens limb 
darkening; the lowest x 2 among them is 41. The Xmin over the entire space of / band models is still 
moderately high at 14. In general, it is the I band that is causing most of the discrepancy between 
models and data. The atmospheric models that have the lowest x 2 form a track in parameter space 
that varies smoothly from solar metallicity dwarfs at T c g = 4750 K to super metal-poor supergiants 
at T cS = 3500 K. 

At this point, we can ask whether the bands are consistent with each other. In a relative sense, 
they are, as the PRMs with the lowest x 2 are always those with T e g = 4500 K, no matter the band. 
We defer discussion as to the possible causes of the disagreement between the ATLAS atmospheric 
models and the microlens data until §4.6, after we have investigated the Next 2 Gen models. 
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4.5. Comparison with Next 2 Gen Models 

We analyze the limb darkening of Next 2 Gen models between T e g = 4000 K and 4600 K in 100 K 
increments, log g = 0.0 and 3.5 in increments of 0.5 and [Fe/H] of —0.25, 0.0 and +0.3. The original 
format of these files is a spectrum between 3500A and 3 /im with a resolution of 0.5A at each of 99 
points in cos 0. Having the full spectra enables us to create limb-darkening profiles in non-standard 
bands. We convolve the spectra with filter functions for all five filters in the microlens data (V, 
Ve, Ie, I and H). We then use a simple linear fit to solve for the (r^,A^) for each filter except H 
band for which we fit for A\ alone. 

The primary difficulty in this procedure is the definition of the edge of a star. The Next 2 Gen 
atmospheres have a steep drop off in intensity whose location in radius varies with surface gravity. 
Sample profiles are shown in Figure 4. This feature cannot be modeled by a linear + square root 
limb-darkening law, and because it contributes almost nothing to the total flux, we decide to remove 
it. This is further warranted because even if we had the formalism in our microlensing code to fit 
this feature, we would not receive any useful information since our sampling is not dense enough at 
the specific part of the caustic exit during which the feature would be visible. We therefore excise 
this feature by removing all points outside of some chosen radius. We then rescale the value of the 
radius at each remaining point by the factor necessary to set the outermost point's radius equal to 
unity. The radius is chosen by finding the point at which the H band (which should suffer the least 
amount of limb darkening), drops steeply off. This is something of a judgment call as individuals 
will pick slightly different truncation radii. This does not pose a problem, however, as tests indicate 
that the Next 2 Gen profile is equally well fit by a two-parameter limb-darkening law out to radii 
outside this steep drop off in flux (but not into the feature we are removing) . In the profiles shown 
in Figure 4, this occurs around r = 0.995, indicating that our cut of r = 0.993 is acceptable. This 
procedure breaks down at low surface gravity. Supergiants have such a small density gradient that 
the surface of last scattering at different wavelengths varies greatly with radius. This would be 
a major concern for us if we did not have additional information telling us that this star was a 
luminosity class III giant. We use a separate truncation radius at each logg, but not for each 
wavelength. Our adopted truncation radius varies between r = 0.88 at logg = 0.0 to r = 0.998 at 
log g = 3.5. Having a single truncation radius for all bands would induce problems at the low log g 
end, but such low surface gravities are already highly disfavored as discussed in §4.3. 

In performing the fit, we sample the profile at the radii corresponding to the observations (see 
Fig. 1), giving equal weight to each point. We evaluate the profile at these radii by interpolating 
among the points given by the Next 2 Gen model. This produces a model profile that is most 
weighted in the regions that are most densely observed in the real data. 

We then perform the same analysis as in §4.4, creating a x 2 surface between the microlens limb- 
darkening parameters and those of the Next 2 Gen models. The only difference is that we now have 
nine limb-darkening parameters instead of five because we are also matching Ve and Ie in addition 
to V, I and H. As in the previous section, no grid points in parameter space for the Next 2 Gen 
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models coincide perfectly with the MPM, and the grid spacing does not perfectly coincide with 
the range of the parameters covered by the PRMs. In this section we will consider T e g- = 4200 K, 
log g = 2.5 and [Fe/H]= +0.3 as the MPM, and we will consider the range T eS = 4100 to 4300 K, 
log g = 2.0 and 2.5, and [Fe/H]= 0.0 to +0.3 as that similar to our PRMs. 

Turning now to the x 2 analysis, the MPM is immediately ruled out: x 2 = 1473. The PRMs 
are also ruled out, as the lowest x 2 among them is 837. When we look over the entire parameter 
space, we find x^m = 227 at T cff = 4600 K, loggr = 2.5 and [Fe/H]= -0.25. From Figure 5(a), it 
appears likely that the true Xmin * s outside of the explored parameter space, at least in T e g. There 
also appears to be a trend towards lower metallicity, so it is possible that the best-fit metallicity 
is also outside our explored parameter space. We will discuss this further in §4.6. In general, x 2 
increases with [Fe/H], and decreases with T c g. Surface gravities between logg = 1.0 and 2.5 are 
preferred, with higher T e Q and lower metallicity selecting for a higher log g. 

4-5.1. Johnson- Cousins V, I and H bands 

At the MPM, the V, I and H bands have x 2 = 93, 237 and 47, respectively. Among the 
PRMs, all three bands prefer T cS = 4300 K and [Fe/H]= 0.0, though the V and H bands slightly 
prefer log g = 2.5 while / band slightly prefers logg = 2.0. The V, I and H bands have x 2 at 
these locations of 32, 172 and 40. The global minima for these bands are: V band) T e g = 4600 K, 
log g = 3.5 and [Fe/H]= 0.0 with x 2 lin = 0.97; I band) T cff = 4600 K, log g = 2.5 and [Fe/H]= -0.25 
with x 2 lin = 60; and H band) T cS = 4600 K, log g = 3.0, [Fe/H]= -0.25 and xL in = 10.3. In sum, 
the all-bands and /-band minima coincide, and the x 2 ' s f° r bands V and H at the all-bands 
minimum is just Ax 2 ~ 2 higher than at their own minima. That is, all three bands have the same 
minimum to within ~ 1 sigma. 

Qualitatively, this shape of the x 2 surface is replicated for each individual band as can be seen 
from Figure 5. The minor differences are that the V band surface has more curvature in T e g, and 
H band has more structure in log g. 

4.5.2. EROS bands V E and I E 

These two bands are in greater disagreement with the atmospheric models than any Johnson- 
Cousin band, except for T c fj = 4600 K and log g > 3.0 (log g > 2.0 for [Fe/H]= —0.25). In most cases 
this disagreement is much greater. The MPM has x 2 = 593 and 503 for Ve and Ie, respectively. 
The best PRM for both bands is T cS = 4300 K, log 5 = 2.0 and [Fe/H]= 0.0 with x 2 = 287 and 302. 
The location of the Xmm f° r the Ve and Ie bands is the same as the all-bands, at T e g = 4600 K, 
log g = 2.5 and [Fe/H]= -0.25 with x 2 = 54 and 99. 
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4.6. Possible Systematic Effects 

Logically, there are only four possible sources for the discrepancy between the models and 
the data: 1) problems with the microlensing data, 2) problems with our analysis of the data, 3) 
problems with the atmosphere models, or 4) incorrect comparison of the models and the data. We 
now argue that (1), (2) and (4) are unlikely. 

4-6.1. Individual Observatories 

We test whether the data from an individual observatory drives the combined solution to 
an unsuitable answer. We rerun our fitting routine five times at the (d,q) of the combined-bands 
solution, each time removing a different observatory's data, the exception being that we always keep 
both H band data sets. The removal of the SAAO, Perth or EROS data sets do not appreciably 
change the limb-darkening curves. In the V band, the removal of the Canopus data shifts the 
limb-darkening parameters by approximately 1 a, and the removal of the YALO data shifts them 
by about 3 a. The reverse is true in the / band, removing the Canopus data provokes a 3 a change 
while removing the YALO data only induces a la change. This test shows that any systematic 
effects in the data themselves are either present in data sets across all observatories, or are so mild 
that they do not affect the combined solution at the level of the difference between microlensing 
measurements and the atmospheric models. 

4.6.2. PLANET vs. EROS Datasets 

We also analyze the solutions found by the PLANET data and EROS data separately. We 
expand the analysis to include searching for a solution over d and q. The PLANET-only solution 
is located at the same (d,q) as the combined solution. Moreover, as discussed in §4.6.1, removing 
the EROS data does not appreciably change the limb-darkening profiles found by the remaining 
(i.e. PLANET) data. However, we find a very different result for the EROS-only solution. This is 
located at a (d,q) of (1.94,0.77) that is (0.005,0.02) away from the combined solution. The EROS 
solution has a x 2 75 lower than the x 2 of the EROS bands at the combined solution. The derived 
stellar brightness profile is flat across the inner half of the star, then drops dramatically towards 
the limb. The ratio of intensities of the center to the limb are similar to that of the V band, but 
the shape of the EROS profiles are very different. The EROS-only profiles are a much better match 
to the Next 2 Gen models than the EROS profiles derived at the all-bands geometry. The x 2 °f the 
MPM drops from 593 and 503 to 174 and 69 for the Ve and Ie bands, respectively. However, the 
profiles from the EROS-only solution are still not actually consistent with any of the Next 2 Gen 
models, since Xmm = 87 and 32 in the two EROS bands. 

Such a major inconsistency is a potentially serious problem. How can the previously described 
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disagreement between atmospheric models and microlensing be trusted if the microlensing can 
produce such different fits? We argue that this problem can be resolved in the following manner: 
1) we identify the feature that has the most diagnostic power with regards to the limb darkening: 
the caustic exit, 2) we show how the EROS datasets do not well determine this feature, though 
the combination of EROS and SAAO H band datasets do, and 3) we investigate whether the large 
formal difference in x 2 between the EROS and all-bands geometries represent a failing of the model 
or of the data, and what the consequences of that failure are. 

We first examine the region from which we receive the most information about the limb dark- 
ening: near the caustic exit. This region is shown in Figure 6. An accurate estimate of the 
caustic-exit time is essential for determining the amount of darkening on the extreme limb. This 
can be illustrated by thinking of data points just outside of the true caustic exit. If the caustic 
exit were thought to occur later than it actually does, these points would be thought to be inside 
of the caustic. Their faintness would therefore imply that the part of the star undergoing the 
strongest differential magnification (the limb) had very low surface brightness. On the other hand, 
if the caustic exit is recognized to occur before these points, their faintness is properly attributed 
to the fact that there are no additional images of the source, i.e. the source lies entirely outside 
the caustic. 

This is exactly the issue with respect to the disagreement between the PLANET-only and 
EROS-only geometries. As Figure 6 shows, the EROS-only geometry places the caustic exit at a 
later time, which implies greater limb darkening. We therefore investigate how well the caustic-exit 
time is determined by the all-bands data set and what the source of this discrepancy is. 

First we note that the time of the caustic exit is essentially determined from the combination 
of the SAAO H band data and the EROS data: the H band data show an approximately linear 
fall toward the caustic exit, and this fall must break very close to the best-fit caustic exit if the 
magnification curve is to remain continuous and still pass through the EROS Ie points. See 
Figure 6. Thus, the caustic-exit time can be specified virtually without reference to any model. 

Even if the EROS points are eliminated from the fit, the best-fit light curve based on PLANET- 
only data still passes through these EROS points and intersects the linearly-falling SAAO H band 
lightcurve at almost exactly the same caustic exit. Because the PLANET points that fix the post- 
exit magnification start up about 10 hours after the exit when the magnification has already started 
to rise, this determination of the caustic exit is somewhat model-dependent. However, the model 
dependence is quite weak. Hence, we have two independent and robust lines of evidence fixing the 
caustic-exit time, and for this reason we have high confidence in the result. 

Nevertheless, it remains somewhat puzzling why the EROS-only solution prefers a later time. 
From Figure 6, it is clear that the EROS data near the caustic exit do not themselves strongly 
prefer one solution over the other. Hence, this discrepancy must be rooted in other parts of the 
light curve: either the EROS data have systematic errors elsewhere in the light curve, or the model 
does not exactly reproduce the true lightcurve. 
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To distinguish between these possibilities, we first exclude the H band data and refit the 
lightcurve. The result is shown by the blue line in Figure 6, which is between the EROS-only 
solution and the PLANET-only solution. Clearly the pressure toward a late caustic exit is not 
coming from the EROS data alone. To verify this, we eliminate the H band, the remaining SAAO 
bands, and the EROS data. The resulting exit, which lies almost on top of the blue line of the 
minus H band geometry, also lies half way between the PLANET-only and EROS-only values. 
Since the problem is not restricted to one data set, we conclude that the model must imperfectly 
predict the data elsewhere in the lightcurve. We have attempted to isolate this discrepancy using 
various techniques, but have not succeeded because the effect is extremely small and only manifests 
itself when data far from the caustic exit are used to predict the caustic-exit time. As with any 
such extrapolation, and caustic-exit predictions in particular (Albrow et al. 1999), small errors can 
be vastly magnified when predicting distant effects. 

We must also determine how much the fit to the caustic region (and so the limb-darkening 
measurement) is being distorted on account of data far away. To do so, we decrease the error by a 
factor of ten on the group of H band points just before the caustic exit. This should increase the 
relative importance of this region to the overall fit and be able to tell us to what extent the light 
curve near the caustic exit is being influenced by data far away. The caustic exit shifts slightly to 
an earlier time, as we should expect given that data far from the caustic tends to shift it to later 
times. This shift, however, does not produce substantial changes in the limb-darkening curves, 
shifting them by 1-1.5 a on average. We conclude that the caustic exit is very well determined, and 
the discrepancies related to it do not significantly affect the limb-darkening determinations. 



4-6.3. Fitting Routine 

It is unlikely that our fitting routine is the source of the conflict. We fit all the data simul- 
taneously, so one would expect any systematic effects would have to be present in all bands. The 
limb-darkening curves in these five bands are all internally consistent with each other. Without 
outside information, the microlensing routine found that the amount of limb darkening increases 
with decreasing wavelength, starting from a very flat profile in H band and progressing through 
/, Ie, and Ve to a relatively steep V band profile (see Fig. 1). Even the relative amounts of limb 
darkening are roughly correct. The mean wavelength of the Ie band is 14% of the way between 
that of the / and V bands (Afonso et al. 2000), and indeed, the Ie band limb-darkening profile is 
very similar to the / band profile. The mean wavelength of the Ve band is 27% of the way between 
V and / (Afonso et al. 2000), and the Ve profile is more like the V band profile than the / band 
profile, as expected the degree of similarity between the Ve and V profiles does not match that 
between the Ie and /. 
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4-6.4- Atmospheric Models 

Figure 7 compares the / band limb darkening derived from microlensing to a suite of ATLAS 
atmospheric models all at the same logg = 2.5 and [M/H]= +0.3, but with different temperatures. 
A striking feature of this figure is that all of the atmospheric models go through a single point, one 
which the microlensing model does not go through. Note that the curve for the Sun (the only star 
in the sky with better measured limb darkening than EROS BLG-2000-5) also passes through this 
point. That is, the models all carry a common feature that is also present in the Sun, but does not 
exist in the giant EROS BLG-2000-5. This exact feature is also present in the Next 2 Gen models 
shown in Figure 8. We conjecture that these models, which are capable of producing fits to the 
Sun, incorporate physics that is only applicable in the solar regime. This common feature is not 
apparent when the limb-darkening curves are normalized in the usual (c\,d\) formalism. It appears 
only when limb-darkening curves are plotted to conserve total flux as in equation (1). In any event, 
this feature gives an easily identifiable handle on the difference between the limb darkening of the 
models and the limb darkening observed via microlensing. 

4-6.5. ATLAS-NexfGen Comparison 

The x 2 surfaces for the ATLAS and Next 2 Gen models are shown in Figures 3 and 5. As 
there are three input parameters, T e g, log g and metallicity for each model, we hold the least 
important constant, and show x 2 versus the other two. This results in two fairly similarly looking 
surfaces. In both cases, T e g is the most important parameter in determining \ 2 with the metallicity 
and surface gravity playing far more subordinate roles. Additionally, while the true Xmin f° r the 
Next 2 Gen models is beyond the range in parameter space we investigated, the flattening of the \ 2 
surface indicates that it is not very far away. Both the ATLAS and the Next 2 Gen models are most 
consistent with a star of T e Q ~ 4700 K, logg = 2.5-3.0 and metallicity around —0.3. 

4-6.6. Model Atmosphere- Microlens Profile Comparison 

As we have just shown in §4.6.5, the atmospheric models both seem to prefer a Kl star. It is 
therefore important to ask if our observations are biased and we are expecting the wrong type of 
star. In principle, differential extinction across the field could redden the star more than the clump 
stars against which it is calibrated (see Fig. 10 of An et al. (2002)), thus affecting our photometric 
estimate of its intrinsic color. However, the spectroscopic analysis yields a source temperature 
similar to the one we find photometrically. 

One should note that reddening has the effect of shifting the observed bands. We find, however, 
that the measured extinction for EROS BLG-2000-5 [E(V-I)=1.3] shifts the mean wavelength of the 
/ bandpass by only +100A. This shift is about 28% of the difference between the mean wavelengths 
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of I and Ie and only 4% of the difference between I and V. Comparing Figures 7 and 1, one sees 
that it cannot account for differences between the observations and the atmospheric models. This 
is partly because the magnitude of the effect is too small, and partly because it causes a shift in 
the wrong direction. 

5. Conclusion 

The observational limb darkening found from microlensing formally disagrees with the limb 
darkening derived from atmospheric models by many sigma. We have argued that this difference 
is unlikely to be the result of the observations, but is more likely due to something related to the 
atmospheric models. It is possible that these models include physics that is not applicable in all 
surface-gravity regimes. It is a testament to the theoretical models that they approximate reality 
in several bands without any previous physical data. We hope that now that giant stars have a 
calibration point, much like dwarfs have from the sun and supergiants have from interferometry, 
stellar models can continue to improve in all stellar regimes. 

We would like to thank Joachim Wambsganss for his comments on this manuscript. DF and 
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Table 1. Data Sets 



Band 


Observatory 


N points 


a scaling 


v 


SAAO 


177 


1.47 




Canopus 


154 


1.96 




YALO 


233 


1.25 


V E 


EROS 


830 


1.59 


Ie 


EROS 


904 


1.63 


I 


SAAO 


404 


1.79 




Canopus 


311 


2.57 




YALO 


424 


1.57 




Perth 


148 


2.89 


H 


SAAO 


549 


2.69 




YALO 


659 


1.95 



Table 2. Model Parameters for EROS BLG-2000-5 





Joint Solution 


An et al. (2002) 


d 


1.935 


1.928 ± 0.004 


q 


0.75 


0.7485 ± 0.0066 


a' 


73? 669054 


74?18±0?41 


u c 


-5.110 x 10~ 3 


(-5.12 ±0.03) x 10- 3 


t> 

L E 


100.030 days 


99.8 ± 1.5 days 


tc 


1736.941667 a 


1736.944 a ±0.005 


P* 


4.792 x 10~ 3 


(4.80 ± 0.04) x 10" 3 


|| 


-0.166358 


-0.165 ±0.042 


7TE,± 


0.195334 


0.222 ±0.031 


d 


0.283 yr- 1 


0.203 ±0.016 yr" 1 




0.0066 rad yr" 1 


0.006 ± 0.076 rad yr" 1 



a Heliocentric Julian Date - 2450000. 
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Table 3. Limb-Darkening Parameters, Errors and Correlation Coefficients 





V band 


Ve band 


Ie band 


/ band 


H band 


r 


0.7077 


0.4813 


0.6494 


0.5173 




or 


0.1021 


0.0855 


0.0656 


0.0726 




A 


0.0456 


0.2880 


-0.3305 


-0.1336 


0.4782 




0.1871 


0.1621 


0.1264 


0.1360 


0.0154 


Cij 


-0.9966 


-0.9979 


-0.9978 


-0.9974 
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0.2 0.4 0.6 0.8 1 

Fractional Radius [ = sin((9)] 

Fig. 1. — The derived brightness profiles for the source star of EROS BLG-2000-5 in all five bands, 
V, Ve, Ie, I and H. The shaded regions around the curves are the 3 a error envelopes. Also shown 
are points indicating when data were taken at each observatory, expressed in distance from the 
caustic to the center of the source in units of the source radius. At this point, the magnification 
pattern over the source profile is discontinuous, which gives us precise information about the profile 
at that point. For example, since the I band has almost continuous coverage across the entire source 
star, we can be confident that the brightness profile well represents reality. The V and H bands 
have a gap in coverage around sin# = 0.8 and thus we can be less confident of the profile there. 
The EROS data, produced by a single observatory, has large gaps, but are still able to make the 
two parameter fit on account of their very dense coverage when the observatory was active. 
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Fig. 2. — The / band x 2 surface plotted as (r, A, \ 2 ) an d viewed along the minor axis of the 
error ellipse. Each point represents a slightly varied geometry and its associated \ 2 an d limb- 
darkening parameters. The filled triangles are those points we use to create the covariance matrix. 
The parabola is the slice through the geometric-error paraboloid for this particular viewing angle. 
Numerical noise will lift points in \ 2 at a particular (T,A), and is responsible for the larger \ 2 
among those points that we do not use. 



All-bands (a) 



V band (b) 




Fig. 3. — Shown are the x 2 surfaces of the microlens- ATLAS comparison at log g = 2.5 for (a) the 
V + I +H bands combined, (b) the V band, (c) the I band and (d) the H band. The value of 
X 2 is most dependent on T e g, and less dependent on [M/H]. The dependence on log g (not shown) 
is even weaker still. The surface of every band has a similar shape, a "valley" that runs through 
metallicity with almost constant T e g. These structures mostly overlap, though the I band's is 
shifted to slightly higher T e g. 
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Fig. 4. — Sample limb-darkening profiles from the Next 2 Gen model corresponding to the MPM 
(T eff = 4200 K, log g = 2.5, [Fe/H]= +0.3). Shown are the V, I E and H bands. Each circle is a 
radial point in the Next 2 Gen output. Also shown is the cut in radius that we impose at r = 0.993. 
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All-bands (a) V band (b) 




Fig. 5. — Shown are the x 2 surfaces of the microlens-Next 2 Gen comparison at [Fe/H]= +0.3 for 
(a) all the bands, (b) the V band, (c) the Ve band, (d) the Ie band, (e) the / band and (f) the 
H band. All the surfaces share the same general shape: monotonic dependence on T c g and [Fe/H], 
with a favored value of log g. With the exception of very low log g (0.0-0.5), \ 2 varies most in T e g, 
with variations over [Fe/H] and logg being much lower. 
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1733.5 1734 1734.5 

HJD - 2450000 

Fig. 6. — The magnification for EROS BLG-2000-5 in the H and Ie bands. The H band has been 
shifted by +1 in magnification to separate it from the Ie band. Also shown are the error bars for 
the SAAO H band and EROS Ie band data points. The magnification for each data point has 
been reconstructed using the observed flux and the source flux, blend flux and seeing correction, the 
last three of which are derived from the all-bands solution. These three quantities are very stable 
and do not appreciably change between microlens models. The black lines show the predicted light 
curve derived from the geometry at the EROS-only solution, the blue line shows the predicted light 
curve from the geometry at the solution containing all data sets except the H band, and the red 
lines show the predicted light curve for the geometry at the all-bands solution. The inset shows an 
expanded view of the Ie caustic-exit region. Note that the Ie band points could support either 
prediction, while the H band points strongly favor the all-bands solution. 
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Fig. 7. — The brightness profile in the I band. The black line and shading are EROS BLG- 
2000-5 and its associated 3cr error envelope. The red lines give ATLAS models (Claret 2000) at 
log g = 2.5 and [M/H]= +0.3 for T eff between 3000 and 10000 K, in 1000 K increments. The blue 
line corresponds to the best-fit ATLAS model at T eff = 4500 K, log g = 3.0 and [M/H]= -0.3. Also 
shown is the I band curve of the Sun in green. As can be seen, there is a point that all of the 
model curves share with the Sun at approximately (r,S7)=(0.74,1.02). This fixed point which is 
not shared by the microlensing-determined profile is not apparent in the (c\,d\) formalism. 
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0.2 0.4 0.6 0.8 1 

Fractional Radius [ = sin(0)] 

Fig. 8. — The brightness profile in the I band. The black line and shading are EROS BLG-2000-5 
and its associated 3 a error envelope. The red lines give Next 2 Gen models at logg = 2.5 and 
[M/H]= +0.3 for T e ff between 4100 and 4600 K, in 100 K increments. The blue line corresponds to 
the best-fit Next 2 Gen model at T eff = 4600 K, log g = 2.5 and [M/H]= -0.25. Also given is the I 
band curve of the Sun in green. As can be seen, there is a point that all of the model curves share 
with the Sun at approximately (r,5/)=(0. 73,1.02). This fixed point is the same as the one that the 
ATLAS models share (but the microlensing-determined profiles do not) and is not apparent in the 
(c\,d\) formalism. 



